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Phenomenological interpretation of the multi-muon events reported by the CDF 

collaboration 
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Abstract 

We present a phenomenological conjecture of new physics that is suggested by the topology and kinematic properties of the 
multi-muon events recently reported by the CDF collaboration. We show that the salient features of the data can be accounted 
for by postulating the pair production of three new states hi, h-z, and hz with masses in the range of 15, 7.3, and 3.6 GeV/c 2 , 
respectively. The heavier states cascade-decay into the lighter ones, whereas the lightest state decays into a r pair with a 
lifetime of the order of 20 ps. 

PACS numbers: 13.85.-t, 14.65.Fy, 14.60.Fg, 14.80.-j, 12.60.Fr 



The Higgs mechanism provides a scheme for the elec- 
troweak symmetry breaking (EWSB). Electroweak pre- 
cision tests (EWPT) indicate that the standard model 
(SM) Higgs boson is light, m h < 186 GeV/c 2 with 
a central value considerably below the lower bound of 
114 GeV/c 2 from direct searches 0. While it is possi- 
ble that one is misled in interpreting the EWPT data in 
terms of a light Higgs mass, the EWSB sector in the SM 
Lagrangian appears to be the most elusive and the one 
most likely to provide experimental surprises [!, 0, 0, Q • 
In Refs. @, S] , the CDF collaboration presents a set of 
studies of multi-muon events. Reference Q uses multi- 
muon events to measure the correlated bb production 
cross section. That study is extended in Ref. [|[ to ad- 
ditional properties of multi-muon events. Based on the 
properties of these events, this paper explores a possible 
conjecture of new physics that might relate the results of 
Ref. 8] to the EWSB mechanism. 

The data set used in Refs. @, S| has been acquired 
with a dedicated dimuon trigger, and consists of events 
containing two central (\r)\ < 0.7) muons, each with 
transverse momentum px > 3 GeV/c, and with invari- 
ant mass larger than 5 GeV/c 2 . References [7, 8] show 
that, when both trigger (initial) muons arise from par- 
ticles that have decayed within the beam pipe of radius 
1.5 cm, their number and kinematic properties are cor- 
rectly predicted by a SM simulation. According to the 
simulation, approximately 96% of the known sources of 
dimuons, such as Drell-Yan, T, Z°, and heavy flavor pro- 
duction, satisfy this condition. The sum of these contri- 
butions (1131090 events) is for simplicity referred to as 
QCD production Q. Reference [|| also reports the ob- 
servation of 295481 events, referred to as ghost events, 
in which at least one muon originates beyond the beam 
pipe. A large fraction of these events is attributed to 
muons arising from n, K, K s , and hyperon decays as 
well as to secondary interactions in the detector volume. 
However, a small but significant fraction of these events 



has characteristics that cannot be explained by known 
processes in conjunction with the current understanding 
of the CDF II detector, trigger, and event reconstruc- 
tion. Reference Q has investigated the rate and kine- 
matic properties of tracks and muons contained in a 36.8 
cone around the direction of each trigger muon in ghost 
events. The nature of these events is characterized by 
four main features. The impact parameter 1 distribu 
tion of initial muon pairs is markedly different from that 
of QCD events. In 36.8 (cos0 > 0.8) cones around the 
initial muon direction, the rate of additional muons and 
charged tracks is significantly higher than that of QCD 
events. The invariant mass of the initial and additional 
muons distributes differently from that expected from se- 
quential semilcptonic decays of hadrons with heavy fla- 
vor. The distributions of the impact parameter of addi- 
tional muons, and of the distance between the pp collision 
point and secondary vertices reconstructed using pairs of 
tracks contained in a cone, do not correspond to the life- 
time of any known particle. 

Although it is not excluded that the features of 
these events can be later explained in terms of known 
sources Q, this paper presents a phenomenological con- 
jecture of new physics that shows that it is possible to 
account for the various features of these events, and also 
suggests additional observations that are further tested 
with the published data In this conjecture, a fraction 
of the ghost events is attributed to pair production of 
three new states h\, /12, and /13 with masses in the range 
of 15, 7.3, and 3.6 GeV/c 2 , respectively. The heavier 
states cascade-decay into the lighter ones, whereas the 
lightest state decays into a r pair with a lifetime of the 
order of 20 ps 0. 

Figure [TJ reproduced from Ref. § , shows the sign- 
coded multiplicity distribution of additional muons con- 
tained in a 36.8 cone around the direction of an initial 
muon in ghost events. In the plot, an additional muon in- 
creases the multiplicity by 1 when of opposite sign and by 
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10 when of the same charge of the initial muon. Leaving 
aside the case in which no additional muons are found, 
an increase of one unit in the muon multiplicity corre- 
sponds in average to a population decrease of approxi- 
mately a factor of seven. This factor happens to coincide 
with the inverse of the r — ► /x branching fraction (0.174) 
multiplied by the 83% efficiency of the muon detector, 
and suggests that these muons might arise from t de- 
cays with a kinematic acceptance close to unity. We use 
Monte Carlo pseudoexperiments to model the shape of 
the muon multiplicity distribution in Fig. [TJ The pseu- 
doexperiments generate 4 t~ + 4 t + leptons and decay 
them into muons with a 17.4% probability. In the pseu- 
doexperiments, initial and additional muons are identi- 
fied with the measured [|| detector efficiencies of 50% and 
83%, respectively. The pseudoexperiment result is shown 
in Fig. [TJ in which the simulated distribution is normal- 
ized to the integral of the data for multiplicity bins higher 
than 10. The comparison of the muon multiplicity distri- 
bution in ghost events with the toy simulation suggests 
that approximately 13200 events contain 8 r leptons in- 
side a cos 9 > 0.8 cone. One interpretation is that they 
are 8-r decays of objects hi that are relatively light with 
respect to the transverse momentum with which they are 
produced. 

Next, we test the data for an obvious consequence of 
the hi —> 8t conjecture. If the hypothesis is correct, 
one expects that, when the hi transverse momentum is 
very large, the 8r decays produce an average of 9.5 tracks 
with pt > 2 GeV/c in a 36.8 cone. We compare the data 
presented in Ref. Q to simulated events which contain 
hi states produced with large transverse momenta. We 
use the pythia Monte Carlo program [l(| to generate 
fictitious pp — ► H — > hihi events followed by hi — > 8r 
decays |ll|. We use this process for convenience, with- 
out implying that the hi states might be Higgs bosons 
or might arise from Higgs decays. Figure [2] compares 
the average track multiplicity in the data and in some 
of the simulations. In order to reject the QCD contribu- 
tion, Reference [|| considers events that contain at least 
three muons in a 36.8 cone. The figure shows the aver- 
age number of tracks with pr > 2 GeV/c contained in a 
36.8 cone around a primary muon as a function of the to- 
tal transverse momentum of the tracks. The asymptotic 
value of the average track multiplicity in the data adds 
support to the conjecture of a hi — > 8t decay initially 
suggested by the muon multiplicity distribution. 

In the assumption that leptons and tracks contained in 
a 36.8° cone around a primary muon arise from hi — > 8r 
decays, the mass of the hi particle is determined by com- 
paring data and simulations in the following invariants 
proportional to m/^ : (a) the invariant mass of all muons 
in a cone when both cones contain at least two muons; 
(b) the invariant mass of all muons in a cone contain- 
ing at least three muons; (c) the muon invariant mass 
for cones containing exactly three muons; (d) the invari- 
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FIG. 1: Sign-coded multiplicity distribution of additional 
muons found in a cos 9 > 0.8 cone around the direction of 
a primary muon in ghost events. The points are reproduced 
from Ref. 0, in which the QCD contribution has been re- 
moved and the distribution has been corrected for the fake 
muon contribution. An additional muon increases the multi- 
plicity by 1 when it has opposite and by 10 when it has same 
sign charge as the initial muon. The first bin represents cones 
without additional muons. As examples, the third bin indi- 
cates cones with 3 muons with charge (H — ) or ( — h +); 

and the 21st bin indicates cones with 3 muons with charge 
(+ + +) or (— — — ). The solid line is the prediction of the 
toy-simulation of a decay into eight r leptons (see text). 



ant mass of muons and tracks for cones containing three 
or more muons and 5 to 6 tracks. Reference Q shows 
that the request of a large number of muons suppresses 
the QCD contribution. In cases when fewer muons are 
requested, the QCD background is larger and has been 
subtracted in Ref. [8j . These invariant mass distributions 
are shown in Fig. [3] The data are compared to simula- 
tions of the process H — > h\hi with m# = 300 GeV/c 2 , 
and m/jj = 15 and 20 GeV/c 2 , respectively. A mass of 
15 GeV/c 2 , close to the invariant mass of 8 r leptons, 
provides a fair modeling of the data, whereas a mass of 
20 GeV/c 2 appears to be too high. 

We investigate if the data are consistent with the hy- 
pothesis of hi pair production by studying the rate and 
properties of events in which two 36.8 cones contain 
a muon multiplicity larger than that of QCD events. 
In the sample of ghost events isolated by the study in 
Ref. [§], there are 27990 ±761 cones containing two or 
more muons, 4133 ± 263 cones containing three or more 
muons, and 3016 ± 60 events in which both cones con- 
tain two or more muons. Figure 2] plots two-dimensional 
distributions of the invariant mass of all muons and of 
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FIG. 2: Average number of tracks in a 36.8 cone around 
the direction of a primary muon as a function of ~}2pt, the 
transverse momentum of all tracks. Data (•), reproduced 
from Ref. [g], are compared to a H — * hihi simulation with 
m/jj = 15, and (o) run = 115 or (o) ran = 300 GeV/c 2 . 
The transverse momentum distribution of the data is different 
from that of the simulations, but the average number of tracks 
has the same asymptotic value, as indicated by the straight 
line. The QCD prediction (■) is based on the few events 
predicted by the heavy flavor simulation, normalized to the 
number of initial muon pairs in the data and implemented 
with the probability that tracks mimic a muon signal. 



the number of tracks contained in each cone for the 3016 
events. Figure [5] shows that the invariant mass distri- 
bution of all muons contained in the 27990 cones con- 
taining at least two muons is consistent with that of the 
3016 events in which both cones contain at least two 
muons [8j. We compare the data to a simulation of the 
process // — > h\h\ generated with the PYTHIA Monte 
Carlo program, in which the h± pair production is medi- 
ated by a photon exchange. In 8549600 generated events, 
we find 15580 events with a pair of initial muons that pass 
the analysis selection of Ref. [8| . These simulated events 
contain 7997 cones with two or more muons, and there 
are 1044 events in which both cones contain two or more 
muons. The ratio 1044/7997 is in quite good agreement 
with that of the data (3016/27990) [H. Based on these 
comparisons, a fraction of the ghost events seems to be 
more consistent with the conjecture that hi states are 
pair produced. 

The dependence of the hi pair production on the Man- 
delstam variable s is studied by comparing data to dif- 
ferent simulations in the following distributions shown in 
Fig. [SI (a) the invariant mass of all muons for events in 
which both cones contain at least two muons; and (b) 
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FIG. 3: Invariant mass, M, distributions of all muons in a 
36.8 cone when (a) both cones contain at least two muons, 
(b) a cone contains three or more muons, (c) a cone contains 
three muons, and (d) of muons and tracks for cones contain- 
ing 5 to 6 tracks and three or more muons. The fake muon 
contribution has been subtracted. The invariant mass dis- 
tributions predicted by a simulation with m^i = 15 (o) and 
m/ii = 20 GeV/c 2 (histogram) are superimposed to the data 
(•) reproduced from Ref. |{|. 
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FIG. 4: Two-dimensional distributions, reproduced from 
Ref. [8J, of (a) the invariant mass, M, of all muons and (b) 
the total number of tracks contained in a 36.8 cone when 
both cones contain at least two muons. 



the invariant mass of all tracks for events in which both 
cones contain at least two muons. The data distribu- 
tions do not show any evidence of resonant production. 
However, the data distributions fall less rapidly than in 
the // — * hihi simulation, a possible indication that 
the hi pair production is not mediated by a photon or a 
gluon. The // — » hihi distribution shown in Fig. [6] cor- 
responds to a production cross section of 3.4 nb. Events 
in which the hi pairs have larger invariant mass are mod- 
eled with a simulation of the process pp — > H — > hih\ 



3 



4000 
5 3000 



(a) 



V 



2 3 4 5 
M (GeV/c 2 ) 



8 300 

CNJ 

o 

« 200 

> 
LU 

100 





(b) 






t 

* 









2 3 4 5 6 

M (GeV/c 2 ) 



FIG. 5: Distributions, reproduced from Ref. |8j, of the invari- 
ant mass, M, of all muons contained in (a) the 27990 36.8 
cones with two or more muons and (b) the 3016 events in 
which both cones contain two or more muons. QCD and fake 
muon contributions have been removed. 
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FIG. 6: Invariant mass, M, distribution of (left) all muons 
and (right) all tracks with pt > 2 GeV/c for events in which 
both cones contain at least two muons. QCD and fake muon 
contributions have been removed. Simulations of the pro- 
cesses // — * h%h\ (dashed histogram), and H — > h\h\ with 
ran = 150 (o) and 300 GeV/c 2 JA) are superimposed to the 
data (•) reproduced from Ref. 



with m H = 150 and 300 GeV/c 2 . The results of these 
simulations, also shown in Fig. correspond to a pro- 
duction cross section of 50 and 35 pb for ran = 150 and 
300 GeV/c 2 , respectively. Using the acceptances mea- 
sured with these different simulations we again estimate 
that approximately 5% of the 295481 ghost events with 
two initial muons can be explained by pair production of 
h\ particles. 

In the assumption that ghost muons are produced in 
the decays of these states, the fits to the high impact 
parameter tail of additional muons in ghost events Q, 
where the distribution is well modeled by an exponen- 
tial function, can be used to estimate their lifetime to 
be 21.4 ± 0.5 ps. If the hi states decay directly into 8 
t leptons with a lifetime of 20 ps, the impact parame- 
ters of muons due to ghost events would be highly corre- 
lated since the r lifetime is negligible compared to that 
of the hi states. In contrast, the study in Ref. Q finds 
a very weak correlation between the impact parameters 
of muons contained in the same cone. We compare data 
from Ref. § to simulations in order to test the more ele- 
gant conjecture of a three-stage decay hi — * 2 h 2 — » 4 ha, 



where /13 is the particle decaying into r pairs. We 
use simulated samples of the processes // — > hih\ and 
pp — > H — > h\h\ with tuh = 300 GeV/c 2 , in which the 
hi states decay into 8 r leptons through a three-stage 
decay. We attribute in turn a 20 ps lifetime to only one 
of the hi, hi, and /13 states. We measure the correla- 
tion between the impact parameters of muons contained 
in a 36.8 cone for the different cases. The correlation 
factor is Pd p d B — 0.39, 0.15, and 0.05 when the lifetime 
is attributed to the hi, h%, and ft. 3 states, respectively, 
whereas it is measured to be 0.03 in ghost events [8j. 
This indirect method provides evidence for the possible 
existence of two additional states, hi and /13. The latter 
state is long-lived and decays into r pairs. Following the 
assignment of a 15 GeV/c 2 mass to the hi state, one ex- 
pects the hi mass to be in the range 7.1 — 7.5 GeV/c 2 
and the mass ~ 3.6 GeV/c 2 . The observed number 
and properties of the ghost events can accommodate the 
additional pair production of at least one of the hi and 
/13 states. When the 20 ps lifetime is attributed to the 
/13 state, only a few percent of the simulated events se- 
lected as in Ref. [8J| survive the additional request that 
both initial muons originate inside the beam pipe. This 
models closely what happens in the data, as reported in 
Ref. 0. 

Because the production mechanism of these hypotheti- 
cal states is not understood, we cannot use the simulation 
to improve the measurement of the h% lifetime reported in 
Ref. [8j . However, we highlight one difficulty that arises if 
ghost events are due to mixed pair production of hi, hi, 
and /13 states. In Ref. Q, the lifetime of the muon parent 
particle has been estimated by fitting with an exponential 
function the muon impact parameter distributions in the 
range 0.5 — 2.0 cm. We have produced simulated samples 
of h n pair production using the processes // — > h n h n and 
pp — > H — » h n h n with m# = 300 GeV/c 2 and n = 1, 2, 3. 
In these simulations, states heavier than the h 3 particles 
cascade-decay into them and ultimately produce events 
with 4, 8, and 16 r leptons in the final state. We gener- 
ated several samples with /13 lifetimes ranging from 10 to 
40 ps. For simulated events due to /13 or hi pair produc- 
tion, analogous fits to the impact parameter distributions 
return the lifetime value used in the event generation. 
However, in the simulation of hi pair production, some 
of the hi cascade-decay products are not relativistic and 
the decay kinematics has greater complexity. As a con- 
sequence, the fits to the impact parameter distributions 
underestimate the lifetime by approximately 30%. 

The identification of r decays into three hadrons would 
provide additional support to the conjecture of new 
physics used to account for a fraction of the ghost events. 
This decay channel could also be used for a complemen- 
tary measurement of the h% lifetime by using the L xy dis- 
tribution of the three-track secondary vertices, where L xy 
is the distance between the secondary vertex and the pp 
collision point projected onto the transverse momentum 
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FIG. 7: Distributions of the (a) invariant mass and (b) dis- 
tance L xy (see text) of three-track systems in events simulated 
with the process pp — > H — > hi hi with mn = 150 GeV/c 2 . 
The three-track systems are produced by single r decays into 
three hadrons. 
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FIG. 8: Distribution, reproduced from Ref. ||, of the dis- 
tance L xy of fit-constrained vertices of three-track systems 
contained in a 36.8 cone around the direction of an initial 
muon for (a) ghost and (b) QCD events. We select cases in 
which angular cones contain only three tracks. 



of the three-track system. Reference [8j has searched the 
data for such secondary vertices, and has verified the de- 
tector response by using identified K s — > 7r + 7r~ decays. 
That study [H reconstructs three-prong secondary ver- 
tices by using tracks with px > 1.0 GeV/c and \r/\ < 1.1 
in a 36.8 cone around the direction of each initial muon. 
Track systems with total charge of ±1 are constrained 
to arise from a common space point. In this study, we 
use similar criteria to search for hadronic r decays in 
simulated events, and compare our result to the pub- 
lished data [8J. In a sample generated with the process 
pp — > H — > h\h\ with rriH = 150 GeV/c 2 , there are in 
average two r decays into three hadrons per event. Ap- 
proximately 8% of these decays, or 0.16 r hadronic de- 
cays per event, survive these selection criteria. Figure [7| 
shows the invariant mass and L xy distributions of the 
three-track systems that are also identified at generator 
level as r decays into three hadrons. 

The simulation also contains 5.5 three-track combina- 
tions per event that pass the same selection criteria, and 
the signal of the 3-hadron r decays is swamped by the 
combinatorial background. However, this signal is com- 
parable to the combinatorial background in a simulated 
sample of h.3 pair production in which a 36.8 cone con- 
tains sometimes one muon and three tracks from the two 
r decays. In this case, the vertices of the three-track 
systems identify correctly the /13 decay vertex and gen- 
erate an excess of events at positive L xy distances. Ref- 
erence Q presents a subsample of ghost events in which 
a 36.8 cone around the direction of an initial muon con- 
tains only three tracks with pr > 1 GeV/c. Three-track 
systems with total charge of ±1 are constrained to arise 
from a common space point. Figure [8] shows the L xy dis- 
tribution for ghost and QCD events Q. Ghost events 
show an excess at positive L xy that is, however, not 
as large as that of QCD events in which most of the 
muon plus the three-track combinations arise from single 
&-quark decays. Figure [9] compares the invariant mass 
of the three-track systems for positive and negative L xy 
values Q. In the case of ghost events, the invariant mass 
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FIG. 9: Distributions, reproduced from Ref. [8J, of the invari- 
ant mass, M, of three-track systems in (a) ghost and (b) QCD 
events. Systems with distance L xy > 0.04 cm (•) are com- 
pared to those with L xy < —0.04 cm (o). Ghost events with 
positive L xy exhibit an excess of events with the expected 
shape of r into three hadron decays. 



of three-track systems with positive L xy exhibits an ex- 
cess of events consistent with the shape of r decays into 
three hadrons shown in Fig. [3 

In conclusion, we suggest one possible phenomenolog- 
ical interpretation of the multi-muon events recently re- 
ported by the CDF collaboration. As shown by the 
comparisons between the published data and simulations 
based on our conjecture, the most interesting features 
of these events can be accounted for by postulating the 
pair production of three new states hi, h 2 , and with 
masses in the range of 15, 7.3, and 3.6 GeV/c 2 , respec- 
tively. The heavier states cascade-decay into the lighter 
ones, whereas the lightest state decays into a r pair with 
a lifetime of the order of 20 ps. The mechanism that 
produces h\ pairs is completely obscure. It does not ap- 
pear to be resonant nor mediated by a photon or gluon 
exchange. The observed pair production cross section 
(~ 100 nb) is a few orders of magnitude larger than what 
is predicted if the h n states belonged to the Higgs sec- 
tor HQ. 
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